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Abstract. We present measurements of the growth rates of the principal facet surfaces of ice 
from water vapor as a function of supersaturation over the temperature range —2 > T > —40 C. Our 
data are well described by a dislocation-free layer-nucleation model, parameterized by the attachment 
coefficient as a function of supersaturation a(a) — Aexp(— ao/cr). The measured parameters A(T) 
and <Jq(T) for the basal and prism facets exhibit a complex behavior that likely originates from 
structural changes in the ice surface with temperature, in particular the onset and development of 
surface melting for T > —15 C. From ao(T) we extract the terrace step energy j3(T) as a function 
of temperature for both facet surfaces. As a basic property of the equilibrium ice surface, the 
step energy j3(T) may be amenable to calculation using molecular dynamics simulations, potentially 
yielding new insights into the enigmatic surface structure of ice near the triple point. 

1 Introduction 

Surface melting occurs when the near-surface atomic or molecular layers of a crystalline solid are not 
as tightly bound as the deeper layers, causing the near-surface layers to lose their ordered structure 
(for a review see p]). The result is an amorphous "premelted" layer, also called a quasi- liquid layer, 
that exists in equilibrium at the solid surface. The structure of the quasi-liquid layer is strongly 
temperature dependent, and its thickness typically diverges as the melting point is approached. 
Surface melting is a common phenomenon in metals and other simple crystalline materials, and it 
has been especially well studied in ice [2j [3j 0] . 

In general the effects of surface melting on crystal growth have been little explored, although 
clearly the structural changes associated with surface melting can have a profound effect on surface 
molecular dynamical processes [5J |B]. Since our theoretical understanding of surface melting is 
relatively poor, we sought to examine how changes in surface structure with temperature affect 
crystal growth behavior. 

It has long been long suspected that surface melting plays an important role in the growth 
dynamics of ice crystals from water vapor |8l |9] . Although ice is a monomolecular crystal with 
a simple hexagonal structure under normal atmospheric conditions, ice crystals forming from water 
vapor exhibit an exceedingly rich spectrum of plate-like and columnar morphologies as a function of 
temperature and supersaturation over the temperature range > T > —40 C [9l [101 [TTJ . Since the 
premelted layer in ice develops over this same temperature range [2j [3l [4] , the prevailing thinking 
holds that the temperature-dependent effects of surface melting on ice crystal growth are responsible 
for the observed morphological complexities, together with instabilities arising from diffusion-limited 
growth and other effects [5]. To date, however, this long-held hypothesis has remained largely 
untested by quantitative experimental data. 
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Figure 1: (a) An idealized schematic of the inner sub-chamber of our experimental set-up. The top 
surface, covered with ice crystals, acts as a water vapor reservoir at temperature Tm for a growing 
test crystal resting on the substrate at temperature T su b s t- When Tm > T su b s t, growth rates are 
determined by measuring the thickness of the crystal (i.e. the distance between the substrate and 
the parallel top facet) as a function ot time using optical interferometry, and by measuring the 
other crystal dimensions using optical microscopy viewing from below the substrate, (b) A rough 
schematic of the actual test chamber, located inside a larger vacuum chamber. The substrate can be 
rotated to bring test crystals into position under the ice reservoir. Additional experimental details 
are provided in [12] . 



As part of our investigation, we made precise measurements of the growth rates of small faceted 
ice crystals from water vapor under carefully controlled conditions, in order to better quantify and 
parameterize the intrinsic ice growth behavior. To this end we measured growth rates of the basal 
(0001) and prism (1100) ice surfaces as a function of water vapor supersaturation a and temperature 
T over the temperature range — 2 > T > —40 C, thus covering the onset and development of surface 
melting. Our measurements were made at low background pressure to reduce the effects of particle 
diffusion through the surrounding gas, so the growth was mainly limited by surface attachment 
kinetics. 



2 Ice Crystal Growth Measurements 

The goal of our ice growth experiments was to examine the growth of individual ice crystal facets in 
a carefully controlled environment, and an idealized schematic diagram of our experimental set-up 
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is shown in Figure [T^. The top surface of the experimental chamber consists of a thermal conductor 
at a uniform temperature Tir, and a layer of ice crystals on its lower surface serves as a source of 
water vapor. At the beginning of each measurement, a single test crystal is placed near the center 
of the bottom substrate surface held at temperature T su b st . The ice reservoir and the substrate 
are separated by thermally insulating walls with a vertical spacing of 1.0 mm. The temperature 
difference AT = Tm — T su b s t determines the effective water vapor supersaturation seen by the test 
crystal. 

During a typical experimental run, we continuously nucleated ice crystals in a much larger outer 
chamber containing ordinary laboratory air (see |12j). where they grew while slowly falling to the 
bottom of the chamber. Typically > 10 7 of these micron-scale crystals were growing within the 
outer chamber at any given time, and the fall times were approximately 3-5 minutes. This cloud of 
slowly growing crystals served as the source of seed crystals for our growth measurements. 

To select a test crystal, an inlet valve on the top of a smaller inner chamber was opened, and air 
carrying suspended ice crystals was drawn from the outer chamber through the inner chamber. The 
operator rotated the substrate while observing the test region under the ice reservoir (see Figure [TJ, 
thus examining crystals that randomly landed on the substrate. When a suitable test crystal was 
identified, the inlet valve was closed and the pressure in the inner chamber was reduced to < 30 
mbar. 

Once the pressure was stable, the operator first adjusted Tir and observed the test crystal 
growing or evaporating slightly in order to determine the a — point, which typically took a few 
minutes. After this, AT was slowly increased to grow the test crystal. The thickness of the crystal 
- defined as the distance between the substrate and the parallel top facet — was determined using 
optical interferometry, while optical imaging was used to record the crystal size and morphology in 
the substrate plane. These data, along with the temperature difference AT, were all recorded as the 
crystal grew. After a few minutes, when the overall crystal size exceeded ~ 100 /im, the crystal was 
evaporated away along with any other crystals on the substrate, and another crystal was selected. 

One important aspect of our experiment was the accuracy of the water vapor supersaturation 
<r seen by each test crystal. Our typical precision in setting AT was ±0.003 C, and drifts in AT 
occurring during a typical measurement were comparable. This overall temperature uncertainty 
corresponded to an uncertainty in the a = point for each crystal of Act w ±0.03 percent. If extra 
care was taken to determine the a — point and stabilize the temperature at the initial stage of a 
growth run, then this uncertainty could be reduced perhaps a factor of two further. 

We checked our calculated supersaturation, based on the measured AT, by nucleating water 
droplets on the substrate (in the absence of any test crystals) and measuring AT at which the 
droplets were neither evaporating nor growing. From this we extracted the supersaturation of water 
relative to ice <J wa ter as a function of temperature, shown in Figure [2] Accurately determining the 
droplet stability point became difficult at the higher temperatures, and this difficultly likely explains 
the systematic trend away from <J wate r(T) seen in Figure O The fact that our measurements, using 
no free parameters, were in good agreement with the accepted cr water (T) suggests that the water 
vapor supersaturation was quite well known in these experiments. 

Only crystal facets parallel to, and thus not in contact with, the substrate were used to determine 
the intrinsic crystal growth parameters. Facets that intersected the substrate often grew at somewhat 
higher rates, especially at low a, owing to substrate interactions that reduced the normal nucleation 
barriers on these facets [15] . 

Broad-band interferometry was used when measuring the basal facets, as described in [121 113j . 
allowing an absolute measurement of the crystal thickness. This technique worked well for thin 
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Figure 2: Measurements of the supersaturation of liquid water with respect to ice in our test chamber, 
determined by nucleating and observing water droplets on the substrate. The line gives the accepted 

CT W ater(T) from [TJ] . 

plate-like crystals, but was less effective when measuring the growth of the prism facets, owing to 
the greater distance between the top facet and the substrate. Thus all the prism facet measurements 
were taken using laser interferometry as described in [T9] . 

For example, Figure [3] shows two images of a single ice prism taken at different times. The central 
laser spot oscillated between dark and bright as the crystal thickness increased, resulting from the 
interference of reflections from the sapphire/ice and ice/vacuum interfaces. These brightness changes 
were used to measure the growth velocity of the top prism facet. For these data the initial crystal 
morphology was assumed to be a simple hexagonal prism, so imaging of the crystal yielded an 
estimate of the initial crystal thickness. 

While our crystal size measurements were straightforward in most cases [12] . we did run into 
some initial problems using laser interferometry on smaller crystals. In addition to the central 
"bulls-eye" fringes seen in Figure El one can also see a number of smaller laser fringes that we 
have come to calling "false" fringes, resulting from multiple reflections within the crystal. The false 
fringes tend to dominate the laser interference pattern when the top facet is small, and the resulting 
bright/dark transitions are not always obvious to interpret. We came to realize that in an earlier 
version of this experiment [19] . we sometimes inadvertently used the false fringes when measuring 
prism facet growth, thus leading to erroneous results. The basal growth measurements in |19] were 
not compromised by this problem, since the basal facet surfaces were large enough to produce clean 
laser fringes. 

Some other important aspects of the experiment included: 1) Our test crystals were small, 
typically < 100 /jm in overall size. Thin plates were typically < 5 /jm thick. The background air 
pressure in the test chamber was typically < 30 mbar. Using small crystals and low background 
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Figure 3: Two images of a single ice prism taken at different times. Interference between reflections 
from the substrate/ice and ice/vacuum interfaces produced the central laser spot seen in these 
images. The brightness oscillated between dark (top image) and bright (lower image) as the crystal 
grew. The c-axis length of this prism is 62 ^m. 

pressures was important to reduce the effects of particle diffusion, as described further below; 2) We 
used only crystals with simple morphologies and well formed facets, and each crystal was discarded 
after growth. Evaporating and regrowing crystals was found to result in generally lower quality data; 
3) Low outgassing materials were used to construct our experimental chamber, and the thermal 
control hardware (including thermoelectric modules, thermally conducting grease, and associated 
wiring) was all mounted outside the vacuum envelope. The chamber was also baked between each 
run to remove any volatile chemical residues. We believe that chemical impurities were responsible 
for some discrepancies in our earlier results [19J . and these problems have been remedied for the 
current experiment; 4) Our test crystals were freshly made in a clean environment and transported 
within minutes to our test chamber with minimal processing, as described above. This also helped 
minimize any buildup of chemical impurities on the ice surfaces. With all these precautions, we 
believe that we have adequately reduced many systematic errors that appear to have been present 
in previous ice growth experiments [9l 116) . 

The goal of our measurements was to determine the intrinsic growth rates of the principal 
facets of ice, which we write in terms of the surface attachment coefficient oti ns t r insic using v = 
o-intrinsicVkinO ' sur / j where v is the growth velocity normal to the surface, Vfcm(T) is a kinetic velocity 
derived from statistical mechanics [H] , and (J surf ace is the water vapor supersaturation immediately 
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above the growing ice surface. The maximum allowed growth velocity has aintrinsic = 1 (albeit with 
some caveats; see [9]), while faceted surfaces generally have OLintrinsic < 1- 

We also found it useful to define a "measured" attachment coefficient a meas derived entirely from 
experimentally measured quantities using v = ot. m eas'Vkin cr oa-i where Coo is the supersaturation far 
from the crystal. To lowest order, the measured a meas (<7oo) is given by 

/ \ ^intrinsic 

(^intrinsic t^oo ) > a^if f 

as described in [SKU], where ay*// ~ 0.15(i?i/i? e / f)(Pi/ P), Reff is an effective crystal radius, P 
is the background air pressure, R± = 1 /im, and P\ = 1 bar. If the measured crystal growth is 
predominantly kinetics limited, then a sur f ~ and a meas (<Joo) ~ ctintrinsic (o'surf) ■ If the growth 
is mainly limited by particle diffusion through the background gas, however, then a meas (<7oo) ~ 
a d iff- 

Examples showing the growth of the basal facets of two ice crystals are shown in Figure IH where 
the measured growth velocities have been converted to a meos ((Too)- At low background pressures, 
we see that the growth is predominantly limited by attachment kinetics, at least for small crystals at 
low supersaturations. For essentially all our data, we found that the growth was well described by 
a layer nucleation model |17) , and to describe the growth we adopted a simplified parameterization 
of the intrinsic attachment coefficient ttmtrinsicip 'sur f ',T) = Acxp(— ao/<7 sur f) where A and <7o are 
parameters that depend on temperature but not on supersaturation, and u sur f is the supersaturation 
at the crystal surface. 

At high background air pressures and high supersaturations, the growth becomes distorted by 
particle diffusion effects. In FigureHJ the solid lines show fits with (A, ao, adiff) = (1, 2.3, 0.15) and 
(1,2.5,0.0075) for the low-pressure and high-pressure crystals, respectively (with supersaturations 
in percent). The intrinsic attachment coefficient ai n trinsic{o sur f) = exp(— 2.3/er stlr /) is also plotted 
as a dotted line in this figure. Note that we have displayed a m eas(<J oo) and ai n trinsic[<Jsurf) in 
the same plot frame in Figure SI even though these can be significantly different quantities. Our 
measurements directly yield a meas (ooo), wn ile our g° a l is to extract ai n t r insic(o 'sur f) from these 
measurements. 

We explored the hypothesis that (^intrinsic may itself depend on pressure, as the background gas 
could affect the surface molecular dynamics that determines extrinsic- If the gas contains chemically 
active components, then these will adsorb on the ice surface and affect its growth, as is well known 
for the case of chemical impurities in air |18j . 

For pure gases such as nitrogen and oxygen, however, as well as clean air, it is generally assumed 
that aintrinsic does not depend on the gas pressure. Comparisons of growth data taken in air at 
one bar and at 25 mbar, like the examples shown in Figure SJ support this conclusion. While the 
diffusion effects depend strongly on pressure, we found no evidence that aintrinsic was changed by 
the presence of the background gas. We therefore made the implicit assumption in our data that 
measurements of aintrinsic made in clean air are equivalent (after accounting for diffusion effects) 
to the ideal case of ice growth from a gas of pure water vapor. 

Although the diffusion distortions are relatively small at the low background pressures used, 
we nevertheless found it important to correct for these effects, and a detailed description of our 
correction procedure is described in [15]. Essentially, the correction produced an estimate of a sur f 
using the known aoo along with the air pressure and measured growth rates of all the facets as inputs. 
From this, the intrinsic attachment coefficient was then obtained from v — ai n trinsicVkin&surf ■ 
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Figure 4: Measurements of the growth of the basal facets of two ice crystals at -15 C, shown as the 
effective condensation coefficient a meas as a function of supersaturation a^ far from the crystal. 
One crystal (dots) was grown in a background pressure of air at 25 mbar, and the other (squares) 
was grown in a background pressure of one bar. The low-pressure crystal shows mainly kinetics- 
limited growth, while the growth at high pressure is mainly limited by particle diffusion when the 
supersaturation is high. Fit lines are described in the text. In addition, the dotted line shows the 
derived intrinsic attachment coefficient ai n trinsic(o 'sur /) • 

This data correction procedure allowed us to remove the residual diffusion effects and convert 
owas(coo) to otintrinsici^surf)- Figure [5] shows a set of basal growth data taken at T = —12 C, 
both without (top) and with (bottom) the diffusion correction. Note that no adjustable parameters 
were use in this conversion. 

Plotting a as a function of a" 1 as in Figure [SJd shows how the parameters A and ao & re coupled 
in our measurements. The slope in Figure [5b gives ao, while the intercept at a^ 1 = is equal to 
A. Correcting for diffusion effects was essential for determining A accurately, since this parameter 
was necessarily arrived at from an extrapolation of the data to a^ 1 = 0. For all our basal growth 
measurements, and for our prism growth measurements with T < — f C, the data were well described 
using A = f, as shown in the example in Figure [5j In these cases we set A = 1 and fit for ao only, 
which reduced the uncertainty in the ao determinations. 

For prism growth measurements at temperatures T > —10 C, it became clear that A = 1 would 
no longer fit the observations. Figure [5] shows one example of corrected growth data at T = —3 C. 
For these measurements we fit both A and ao to the data. We found that ao was low with these 
data, so coupling between the parameters was not a serious problem, allowing both parameters to 
be measured reliably. 

In addition to a much lower A, Figure [6] also demonstrates small systematic variations between 
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Figure 5: (a) A plot of the measured a mea s ipao) of 21 ice crystals at T = — 12 C. The solid line gives a 
fit a meas = amtrinsicatdif f I ' (^intrinsic + «<«//)> where a intr insic = exp(-0.021/er) and a di f f = 0.07. 
(b) A plot of ctintrinsici.o'surf) for the same crystals after correcting for diffusion effects. The solid 
line gives the fit aintrinsic = exp(— 0.0195/u), and this same function is shown as a dashed line in (a). 
Note that the overall horizontal shift in the corrected points reflects the fact that (J surface < Coo- 

three different data sets acquired during three separate measurement runs. The high and low runs 
(open points) were done while testing different substrate surface treatments, while the central run 
(solid points) was done with a clean sapphire substrate. We suspect that the variations seen in 
Figure |6] resulted mainly from residual chemical vapors in the test chamber, associated with the 
surface treatments, depositing on the top ice surface and affecting its growth. Nearly all our data 
were taken using a substrate that was thoroughly cleaned before each run, with no additional surface 
treatments. 

Residual chemical impurities of this nature are an unavoidable uncertainty in all ice growth 
experiments, since one cannot be certain how clean is clean enough. Furthermore we have seen 
situations were chemical impurities blocked the surface growth and thus reduced aintrinsic, as well 
as circumstances where impurities reduced the nucleation barrier and thus increased aintrinsic- We 
have used quantitative experiments to examine how surface chemistry can affect ice crystal growth 
|18j , and we believe that the current experiment is sufficiently clean that residual chemical effects are 
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Figure 6: Measurements of ai n trinsic(&surf) versus (T~ ur f obtained from 19 crystals growing at a 
temperature of —3 C. The three different sets of symbols are from separate measurement runs done 
on different days. The lines show a — Aexp(— <Jo/a) with parameters (A,<to) = (0.17, 0.00013), 
(0.13,0.0002), and (0.085, 0.00025). 

not a serious systematic problem. The final error bars in our measurements were adjusted to account 
for residual systematic errors of this kind, as estimated empirically from run-to-run variations seen 
throughout all our data. 

For both ice facets, and over the entire temperature range measured, we found that our data 
are well described by the functional form ai n t r insic(o'surf,T) — Aex£>{—uo/<7 sur f), indicative of 
a dislocation-free layer-nucleation model. From measurements of over 200 crystals, we produced 
the final measurements of <Jq(T) and A(T) shown in Figure [3 which is the principal result from 
this experiment. The current set of measurements includes only data points for T > — 20 C; lower- 
temperature points in Figure[7]were taken from |19j . which was a previous version of this experiment. 
For those cases where we assumed A — 1 in our fits, as described above, we assigned error bars to A 
in Figure [7] that give an estimate of the overall experimental uncertainty in this assumed value. The 
ranges of surface supersaturations over which data were collected for these measurements are shown 
in Figure [H 

Roughly 5-10 percent of the crystals sampled grew much more rapidly than the norm, especially at 
low supersaturations, suggesting that the usual nucleation barrier was substantially reduced [TBI 120) . 
We suspect that dislocations or perhaps isolated surface chemical impurities affected the growth of 
these crystals, and they were discarded from our data set before fitting to produce the measurements 
in Figure [7] 

Another set of crystals, again roughly 5-10 percent of those sampled, grow anomalously slowly, 
and these crystals were also discarded from our fits. In essence, we performed "robust" fits to our 
data - first fitting the entire data set, then removing a small number of "outlier" crystals before 
redoing the fits. This step was necessary because the outlier crystals grew substantially differently 
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from the norm and thus adversely affected the fits. 

We have witnessed similar outlier effects in other ice growth experiments, suggesting that this 
is a persistent and somewhat uncontrollable problem. Managing this problem was possible only by 
measuring a large number of crystals to define the overall distribution of growth behaviors, thus 
allowing us to identify and remove highly unusual cases. We believe that the final results in Figure 
[7] accurately represent the growth of chemically clean, dislocation- free ice facet surfaces. 

Because the crystal growth we observed was everywhere well described by a layer-nucleation 
model, the measured supersaturation parameter oq can be used to calculate the terrace step energy 
,9(T) as a function of temperature for both facet surfaces using 

where Q2 is the area of a molecule on the surface. This relation comes from classical 2D nucleation 
theory [17l H9] . A plot of f3(T) from our data is shown in Figure El We note from the scale on the 
right side of Figure [9] that /3(T) is much smaller than /3q = ja « 3.5 x 1CP 11 J/m, the product of 
the surface energy 7 « 0.11 J/m 2 of the ice/ vapor interface and the nominal molecular step height 
a ~ 0.32 nm, which is an upper limit on the step energy |19) . 

3 Discussion and Interpretation 

A comparison of the results from the present experiment with analogous results from earlier ice 
growth experiments is beneficial for advancing the state-of-the-art for these measurements. A careful 
scrutiny of the published results reveals that in many older experiments a number of systematic errors 
were not adequately managed (for a summary, see [HI [16]). Diffusion effects especially distort the 
ice growth behavior, and this fact was not always adequately appreciated in prior experiments. 

Removing diffusion effects to reliably determine oiintrinsic is especially difficult in experiments 
done at pressures near one bar. Precise diffusion modeling would be necessary to extract cxintrinsici 
and we believe that sufficiently accurate modeling techniques have not yet been demonstrated. Fur- 
thermore, our analysis shows that diffusion effects are important even at quite low pressures, and that 
the corrections become larger with larger crystals and at higher supersaturations. By identifying, 
reducing, and modeling this and other systematic effects, we believe that the current measurements 
are a substantial improvement over previous attempts to determine Oiintrinsic(&surf,T). 

Interpreting our results, summarized by the parameterization of cxi n t r i ns i c (a sur f,T) shown in 
Figure U\ presents a significant challenge. Our theoretical understanding of the surface structural 
changes that accompany surface melting is itself rather poor. How these structural changes in turn 
affect the crystal growth dynamics is clearly a complex many-body problem. We offer the following 
conclusions and observations: 

1) Our first conclusion, as stated above, is that the measured attachment coefficients for both 
principal facets are well described by a dislocation-free layer-nucleation model with the simplified 
parameterization cti n trinsic(Psurf ,T) = Aexp(—ao/a sur f). We find it quite remarkable that for both 
facets the growth dynamics of can be summed up so concisely by the functions A{T) and ao(T). 
This is true even going through the transition from essentially no significant surface melting to a 
fully developed quasiliquid layer. Although the equilibrium structure of the ice surface changes 
dramatically over this temperature range, as does the equilibrium vapor pressure, the functional 
form a ~ Aexp(— (To/a) remains unchanged as the ice growth is everywhere described by a layer 



10 




■ 1 1 1 1 I I I I I 1 1 1 1 — h 

- f— HffiffiH-H i 




1 10 

- Temperature (C) 

Figure 7: Measurements of the growth behavior of the basal (solid points) and prism (open 
points) facet surfaces of ice crystals. The intrinsic attachment coefficient was parameterized by 
ccintrin,sic(&surf>T) — A(T) exp(— <tq (T)/<7 sur f), and the plots show the parameters A(T) and uq(T) 
extracted from our data for both principal facets. 

nucleation model. The dominant change with temperature is seen in the nucleation parameter <jq(T), 
accompanied by a relatively modest change in A(T) on the prism facet. 

2) The measured 0o(T) immediately yields the terrace step energy /3(T) from an application of 
classical 2D nucleation theory, as described above. We note that the step energy is a fundamental 
property of the ice surface, in much the same way that the surface energy is a fundamental quantity. 
It is also an equilibrium property, even though it was derived here from the dynamical process of crys- 
tal growth. As an equilibrium, molecular-scale quantity, the step energy /3(T) may be amenable to 
calculation using perturbation techniques or molecular dynamics simulations. Considerable progress 
has been made in investigations of ice surface melting using molecular dynamics simulations |21[ I22) , 
so perhaps step energies can be calculated using similar methods. The observed strong temperature 
dependence in j3(T) may thus yield important insights into how surface melting affects the ice surface 
structure, in particular the interface between the crystalline solid and the quasi-liquid surface layers. 
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Figure 8: The ranges in a sur f over which growth data were collected for the two facets. The ranges 
were set mainly by the minimum and maximum growth velocities we could reliably observe with 
this experiment. 

The step energy can also be used to infer some features of the equilibrium crystal shape, which to 
date has not been reliably measured for ice [23] . 

3) The growth behaviors of the basal and prism facets are remarkably different, especially at 
temperatures above T m — 5 C. In particular, at the highest temperatures we see ao.basai 3> &o, prism, 
in a temperature regime where we would expect to find a thick quasi-liquid layer on the ice surface. 
This is consistent with the growth behavior of ice from liquid water, where the basal surface grows 
much more slowly than the prism surface at low supercoolings. Quantitatively relating the crystal 
growth rates from water vapor and from liquid water may be a tractable theoretical problem, since 
the surface attachment kinetics is likely similar for these two cases. To our knowledge, however, 
the kinetic coefficient for ice growth from liquid water has not yet been determined, as growth 
measurements are usually limited by heat diffusion effects. 

4) For both facets we see an overall trend that j3(T) decreases with increasing T. As the temper- 
ature increases, we suggest that surface restructuring smooths out the terrace edge, thereby lowering 
the step energy. The trend with temperature is particularly simple on the basal facet for T < — 15 
C, where we expect the complicating effects from surface melting are small. These low tempera- 
tures may be best suited for molecular dynamics simulations, so extending the prism facet data to 
similarly low temperature would be beneficial. 

5) The small bump in cro,basai(T) at T » —12 C is a significant and robust feature in our data, 
and we suggest that this feature identifies the onset of significant surface melting (significant with 
regard to its effect on crystal growth dynamics) on the (0001) surface at this temperature. The 
mechanism that might produce this bump, however, is unknown. This onset temperature agrees 
with that measured by [4] , which is a surface melting measurement using an ice surface preparation 
similar to that used in the current experiment. The observed shoulder in (Jo jPr ism{T) near —12 C 
may be a related phenomenon. More generally, we see that the behavior of <jq(T) for both facets is 
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Figure 9: The step energy f3(T) extracted from our measurements of <J (T) using classical nucleation 
theory, for the basal (solid points) and prism (open points) facet surfaces. The scale on the right 
compares j3 with /3q = 7a, where 7 is the surface energy and a is the step height. 

most complex in the temperature range — 15 C < T < — 5 C, and we suggest that this range includes 
the onset and development of surface melting on both facets. At lower temperatures surface melting 
is largely absent, while at higher temperatures the quasi-liquid layer is fully developed. 

6) The change in <jQ i,asai{T) near T = —5.5 C is especially noteworthy, as it reverses the general 
trend with temperature seen in <Jq(T) for both facets. One possibility is that surface melting initially 
results in a more disordered basal surface with a reduced step energy, producing the dip in ao t basai(T) 
near T = — 7 C. At higher temperatures, however, we suggest that the quasi-liquid layer (QLL) 
becomes more fully developed, producing a relatively sharp QLL/ice interface. This sharper interface 
may then result in a more distinct terrace step and a higher step energy. 

7) We note that our measurements of ai n t r insic(o'surf,T) do not lead one immediately to an 
explanation of the well-known morphology diagram describing ice growth from water vapor [9]. In 
fact, the opposite is true. At —15 C, for example, we find otbasai > ciprism at all a, which is at 
odds with the occurrence of very thin plate-like crystals at this temperature. The explanation for 
this and other morphological discrepancies seems to lie in the fact that our measurements give 
o>intrinsic{vsurf,T) only for flat facet surfaces. The edge of a thin plate, however, is not a flat facet 
surface. The molecularly flat prism facet in this case is typically only some hundreds of molecules 
wide (assuming an edge radius of curvature of approximately one micron), and it is the dynamics on 
this surface that determines the edge growth velocity. Earlier one of us put forth the hypothesis of 
structure- dependent attachment kinetics, suggesting that a changes depending on the local crystal 
morphology, and in particular that <jq is reduced on thin edges [24] . We believe that including this 
hypothesis can connect the present measurements with the morphology diagram, although many 
details remain unknown at present. 

In summary, we have measured the ice crystal growth velocities of the principal facets of ice from 
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water vapor for what approaches the ideal case - that of perfect crystalline surfaces growing in near 
equilibrium with pure water vapor. We find that data such as these can be used both as a probe of 
the temperature dependence of surface melting and as a measure of the effects of surface melting on 
crystal growth dynamics. This work was supported in part by the California Institute of Technology 
and the Caltech-Cambridge Exchange (CamSURF) program. 
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